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The invention relates to a computer tomography method in which a bundle of 
rays passes through an object that moves periodically, in particular a heart. The invention 
also relates to a computer tomograph for carrying out the method and to a computer program 
for controlling the computer tomograph. 
5 In known methods of the kind specified in the paragraph above, the 

distribution in space of the absorption or attenuation of the radiation in objects moving 
periodically is reconstructed from measured values that are acquired with a detector unit. 
When this is done, the periodic movement of the object means that the measured values 
convey information from different states of the object, which leads to motion artifacts in the 

10 dataset that is reconstructed. 

To reduce such motion artifacts, what is done in known methods is for 
movement signals dependent on the movement of the object, such as an electrocardiogram 
for example, to be recorded simultaneously during the detection of the transmitted rays. The 
different successive phases of the movement of the object, which are repeated periodically, 

15 can be determined from these movement signals. The measured values that are then used to 
reconstruct the object are solely ones that were acquired when the object was in the same 
phase of its movement. What a phase of movement is in this case is generally a range of 
phases that is repeated in each cycle and in which the times of acquisition of the measured 
values used for reconstruction are situated. If, for example, the phase of movement that has 

20 been selected in the case of a heart is 75% to 85% R-R (R-R = the interval between adjacent 
R peaks in an electrocardiogram), then the measured values that are used for reconstruction 
are solely ones whose times of acquisition were situated in this range of phases in each of the 
individual cycles. 

What is disadvantageous in this case is that, due for example to the varying 
25 rate of heart-beat, the states of the object in the phases of movement in the different cycles 
are not really the same. In this way, in a phase of movement of 75% to 85% R-R in a period 
or cycle, the heart, as an object, assumes states that differ from its states as an object in the 
same phase of movement in a succeeding period, which produces severe motion artifacts. 
This effect is all the more pronounced the more the object moves in the given phase of 
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movement. An attempt is therefore made in known methods of heart imaging to determine, 
with the help of the electrocardiogram, that phase of movement in which the object moves 
least. This however can only be done to an inadequate degree, because the electrocardiogram 
does not reproduce the movement of the object with sufficient accuracy. 
5 It is an object of the present invention to specify a computer tomography 

method, a computer tomograph and a computer program in which the motion artifacts are 
less pronounced. 

This object is achieved in accordance with the invention by a computer 
tomography method having the following steps: 
10 a) generation by a radiation source of a bundle of rays that passes though an 

object moving periodically, 

b) production of a relative movement between the radiation source on the one 

hand and the object on the other hand, which relative movement comprises rotation about an 
axis of rotation, 

1 5 c) acquisition, with a detector unit and during the relative movement, of 

measured values that depend on the intensity in the bundle of rays on the far side of the 
object, 

d) sensing of a movement signal dependent on the movement of the object with a 
movement-sensing means and determination of periodically repeated phases of movement 

20 with the help of the movement signal sensed, 

e) reconstruction of a plurality of intermediate images of a region of the object, 
each intermediate image being reconstructed with measured values that were acquired while 
the object was in a different phase of movement, thus enabling a phase of movement to be 
assigned to each intermediate image, 

25 f) determination of the phase of movement in which there was least movement of 

the object in the region, by determining the intermediate image having the fewest motion 
artifacts in the region, 

g) reconstruction of a computer tomographic image of the region from measured 

values that were acquired while the object was in the phase of movement in which there was 
30 least movement of the object in said region, the reconstruction parameters that are used in 
this case differing from the reconstruction parameters used to reconstruct the intermediate 
images. 

Compared with the known methods mentioned above, what is done in 
accordance with the invention is for that phase of movement first to be determined in which 



WO 2005/055829 



PCT/IB2004/052612 



3 

the object moved least, in the region concerned, during the acquisition. The measured values 
than are then used for the reconstruction are solely ones whose times of acquisition were 
situated in those phases of movement, in which case, as has already been mentioned above, 
the phases of movement are generally ranges of phases. This produces a reduction in motion 
5 artifacts as compared with known methods. 

The term "periodical movement" is not confined to what is meant by a 
periodicity in the exact sense, i.e. it is not confined to movements in the course of which 
states of an object are regularly repeated exactly, i.e. exactly identical states of the object 
occur at exactly equidistant points in time. A periodical movement in the context of the 

10 invention covers in particular movements that are not of mathematical exactness, as is known 
to be the case with organs of the body that move periodically, such as the heart for example. 
What is meant is that there is a progression through similar, substantially identical, states of 
the object at points of time that are generally substantially equidistant. 

The region within the meaning of claim 1 may comprise the whole of the 

1 5 object or only part thereof. 

The reconstruction of the intermediate images with a low spatial resolution as 
claimed in claim 2 leads to a reduction in computing work and cost when the intermediate 
images are being reconstructed. The spatial resolution of the intermediate images has to be 
sufficiently high to enable motion artifacts to be discerned in the intermediate images. No 

20 other requirements are normally laid down for the resolution of the intermediate images. In 
particular, it is not necessary to reconstruct intermediate images of as high a spatial resolution 
as possible such as are required for diagnostic purposes. 

In claim 3 a region of the object that is to be examined (the field of view - 
FOV) is divided into sub-regions. Then, there is determined for each sub-region the 

25 respective phase of movement in which the particular sub-region moved least during the 

acquisition. What are then used for reconstructing a sub-region in the computer tomographic 
image (CT image) to be reconstructed in conclusion are solely measured values that were 
acquired while the object was in a phase of movement in which it moved least in the sub- 
region concerned. Because the object may move differently in different sub-regions, this 

30 results in a further reduction in the motion artifacts. 

Determining the intermediate image having the fewest motion artifacts by 
means of a motion-artifact metric and a motion-artifact value as claimed in claim 4, and in 
particular by means of a mean of gradients of image values of an intermediate image in the 



WO 2005/055829 



PCT/IB2004/052612 



4 

direction of the axis of rotation as claimed in claim 5, leads to a further reduction in the 
motion artifacts. 

When known methods are being used, motion artifacts occur to a particularly 
marked degree at points where regions of the object are reconstructed with measured values 
5 whose times of acquisition, although they were situated in the same phase of movement, were 
situated in different periods. Such regions are referred to as overlap regions. Overlap regions 
occur especially frequently with objects that move fast relative to the data acquisition 
process, such as the human heart, because high temporal resolution is needed to reconstruct 
such objects, which means that ranges of phases of movement that are as narrow as possible 
10 have to be used. Because the measured values used are solely ones whose times of 

acquisition were situated in these ranges of phases of movement then, in order to have a 
sufficient number of measured values available for the reconstruction, measured values from 
as many different periods as possible are used to reconstruct the same region of the object. 
The embodiment claimed in claim 6 allows for this fact by assigning to gradients situated in 
15 these overlap regions a greater weight than to gradient that are not situated therein. This gives 
a further reduction in the motion artifacts. 

A computer tomograph for performing the method according to the invention 
is defined in claim 7. 

Claim 8 defines a computer program for controlling a computer tomograph as 
20 claimed in claim 7. 

These and other aspects of the invention are apparent from and will be 
elucidated with reference to the embodiments described hereinafter. 



25 In the drawings: 

Fig. 1 shows a computer tomograph with which the method according to the 
invention can be performed. 

Fig. 2 is a flow-chart of the method according to the invention. 

Fig. 3 is a schematic representation of an electrocardiogram, a period, and a 
30 plurality of phase points in the period, 

Fig. 4 is a flow chart of a filtered back-projection and 

Fig. 5 is a schematic perspective view of a helical trajectory, a virtual detector 
and a plurality of fans of rays. 
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The computer tomograph shown in Fig. 1 comprises a gantry 1 that is able to 
rotate about an axis of rotation 14 extending parallel to the z direction of the coordinate 
system shown in Fig. 1 . For this purpose, the gantry 1 is driven by a motor 2 at a preferably 
5 constant, settable angular speed. Fastened to the gantry 1 is a radiation source S, preferably 
an X-ray source. This latter is provided with a collimator arrangement 3 which marks out 
from the radiation generated by the radiation source S a conical bundle of rays 4, i.e. a bundle 
of rays that is of a finite extent other than zero both in the z direction and in a direction 
perpendicular thereto (i.e. in a plane perpendicular to the axis of rotation). As an alternative, 

1 0 a fan-shaped beam could also be used. 

The bundle of rays 4 passes through a cylindrical examination region 13 in 
which is situated an object (not shown) that moves periodically. In the present embodiment, 
this object is a beating heart that performs spontaneous movements and that may also, under 
certain circumstances, be moved to and for by respiratory movements of the patient. In other 

15 embodiments, the rays might also pass through other parts of the body that move 
periodically, such as the liver, brain or veins, or though technical objects that move 
periodically. 

Having passed through the examination region 1 3, the bundle of rays 4 
impinges on a detector unit 1 6 fastened to the gantry 1 , which detector unit 1 6 has a detector 

20 surface comprising a plurality of detecting elements that in this embodiment are laid out in 
rows and columns in the form of a matrix. The columns of detecting elements preferably 
extend parallel to the axis of rotation 14. The rows of detecting elements are situated in 
planes perpendicular to the axis of rotation, in the present embodiment along arcs about the 
radiation source S (making the detector unit a focus-centered detector). In other embodiments 

25 however, they may also be differently arranged and may for example describe arcs about the 
axis of rotation 14 or be in a straight line. Each detector element on which the bundle of rays 
4 impinges gives, in each position of the radiation source, a measured value for one ray from 
the bundle of rays 4. If a fan-shaped bundle of rays is used in other embodiments, then the 
detector unit could have only a single row of detectors. 

30 The included angle a™™ of the bundle of rays 4 determines the diameter of the 

object cylinder within which the object being examined is situated when the measured values 
are acquired. The included angle is defined in this case as the angle that a ray that is situated 
at the edge of the bundle of rays 4 in a plane perpendicular to the axis of rotation 14 makes 
with a plane defined by the radiation source S and the axis of rotation 14. The examination 
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region 13, i.e. the object or the patient presentation table, can be moved by means of a motor 
5 parallel to the axis of rotation 14, i.e. to the z axis. As an equivalent to this, it would 
however also be possible for the gantry to be moved in this direction. When the object 
involved is a technical one and not a patient, it may be the object that is turned in the course 
5 of an examination while the radiation source S and the detector unit 16 remain stationary. 

By means of the motors 2 and 5, the radiation source S and the detector unit 16 
are able to follow a trajectory relative to the examination region 13 than extends along the 
surface of an imaginary cylinder. This trajectory may for example be helical in form if both 
motors operate. If on the other hand the motor 5 for feed in the direction of the axis of 

10 rotation 14 remains stopped and the motor 2 causes the gentry to rotate, a circular trajectory 
will be obtained for the radiation source S and the detector unit 1 6 relative to the examination 
region 13. It will be the helical trajectory that is considered in the present embodiment. 

During the acquisition of the measured values, the movement of the heart is 
recorded in a known fashion by means of an electrocardiograph 8. For this purpose, the chest 

15 region of a patient is connected to the electrocardiograph 8 via electrodes (not shown). 

Alternatively, the pulse could also be used as a signal defining the movement of the heart. In 
other embodiments, and particularly with other moving objects, the movement of the object 
can be followed with the help of other movement signals. In this way, with a technical object 
that is moved periodically by a motor, a signal from the motor may be used as a movement 

20 signal. 

In the present embodiment, it is assumed that the patient does not breathe 
during the measurement. The respiratory movements can thus be ignored. Alternatively, the 
respiratory movement could be measured with, for example, a deformable chest belt that was 
connected to a means of measuring respiratory movement. 

25 The measured values acquired by the detector unit 1 6 are fed to a 

reconstructing unit, and particularly a reconstructing computer 10, which is connected to the 
detector unit 16 by for example a wireless data transmission system (not shown). The 
electrocardiogram too is transmitted to the reconstructing unit 10, from the 
electrocardiograph 8. The reconstructing unit 10 reconstructs the distribution of absorption in 

30 the examination region 13 and passes it on, to a monitor 1 1 for example. The two motors 2 
and 5, the reconstructing unit 10, the radiation source S, the electrocardiograph 8, the 
transmission of the measured values from the detector unit 16 to the reconstructing unit 10 
and the transmission of the electrocardiogram from the electrocardiograph 8 to the 
reconstructing unit 10 are controlled by the control unit 7. 
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In other embodiments, the measured values acquired and the 
electrocardiograms measured may be fed for reconstruction purposes first to one or more 
reconstructing computers, which then pass on the reconstructed data to an image-processing 
computer via, for example, a fiber optic cable. 
5 Fig. 2 shows the flow of a measuring and reconstructing method that can be 

performed with the computer tomograph shown in Fig. 1 . 

After initialization in step 101, the gantry rotates at an angular speed that is 
constant in the present embodiment but may also vary, e.g. as a function of time or of the 
position of the radiation source. 

10 In step 102, the examination region, i.e. the object or the patient presentation 

table, is moved parallel to the axis of rotation and the radiation from the radiation source S is 
switched on, thus enabling the detector unit 16 to detect the radiation from a plurality of 
angular positions. As the same time as, or even before, the switching on of the radiation 
source S, the electrocardiograph 8 is activated, which means that an electrocardiogram 2 1 

15 (Fig. 3) is measured at the same time. 

In step 103, phases of movement, which represent ranges of phases in the 
present embodiment, are arranged at different phase points q>\ ...<pn within each period with 
the help of the electrocardiogram 21. For this purpose, one complete period, in which the 
object passes through each phase of movement once, is defined by the interval between two 

20 adjacent R peaks 23 in the electrocardiogram 21 . The phase points cp\ ... q>N are then 

distributed equidistantly for example in each period, e.g. at the points representing 2% R-R, 
4% R-R, . . ., 98% R-R. The expression "x% R-R" designates in the present case a point t = tn 
+ Q.Olx&tR.R, where (r is the time of the first R peak in the period in question and AtR.R is the 
interval of time between the two R peaks in the period in question. The phases of movement 

25 are positioned in each period in such a way that each phase point <p\ ... <p^ is in the center of a 
phase of movement. The width of the phases of movement is determined in the course of step 
104. 

In step 104, an intermediate image of an object region to be reconstructed (a 
field of view - FOV) is reconstructed. What this means is that the measured values that are 
30 used for reconstructing the intermediate image to be assigned to the phase point (p\ are solely 
ones that were acquired while the object was in that phase of movement in the different 
periods that was located at the phase point <p\. The same is true of the intermediate images 
that are assigned to the other phase points. The width of the particular phase of movement is 
adjusted in this case in such a way that a presettable quantity of measured values are 
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available for reconstructing the particular intermediate image. If for example each voxel of 
the FOV is to be irradiated over an angular range of at least 1 80° then, on the basis of the 
known acquisition geometry, of the angular speed and of the electrocardiogram, appropriate 
widths have to be determined for the phases of movement by simple geometrical calculations 
5 or by computer simulations. 

The reconstruction of the intermediate images in step 104 may for example be 
performed by means of a filtered back-projection, which is explained in detail below in 
connection with Fig. 4, or by means of known iterative procedures such as ART (algebraic 
reconstruction technique). Before the reconstruction, each measured value may be weighted 
10 by being multiplied by a weighting factor. This weighting factor may be all the larger, the 
more the corresponding measured value is located in the center of the particular phase of 
movement. 

The reconstruction of the intermediate images is preferably performed with a 
spatial resolution that is as low as possible (e.g. 20 x 20 x 20 cm 3 represented by 32 3 voxels) 
15 but still allows motion artifacts to be detected, the purpose being to minimize computing 
work and costs. 

In step 105, a motion-artifact value is determined for each intermediate image 
by means of a motion-artifact metric, solely from image values from the particular 
intermediate image. 

20 It is known that reconstructed CT images of an object that have pronounced 

motion artifacts have a higher mean gradient of image values in the direction of the axis of 
rotation 14 (z gradient) than CT images of the same object that have fewer motion artifacts. 
The mean of the z gradients of image values can therefore be used as a metric for motion 
artifacts. The smaller this mean is, the fewer motion artifacts there are in the intermediate 

25 image in question and the less did the object move in the corresponding phase of movement. 
A similarity value can therefore be determined by, for example, calculating the z gradient for 
each voxel of the particular intermediate image and finding the means of the z gradients. 

It is also known that a reconstructed CT image of an object that has 
pronounced motion artifacts has image values that vary greatly with changing angular 

30 positions of the radiation source S, whereas in CT images of the same object in which the 

motion artifacts are less pronounced these changes as a function of the angular position of the 
radiation source occur to a less marked degree. A mean gradient in the direction of the angle 
of rotation of the radiation source S can therefore also be used as a motion-artifact metric. A 
motion-artifact value may for example then be determined by calculating a gradient of the 
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image values in the direction of this angle of rotation for each voxel in the particular 
intermediate image and finding the mean of these gradients. The resulting mean is then the 
motion-artifact value for this intermediate image. 

When a motion-artifact value is determined with the help of the z gradient, 
5 each z gradient may be multiplied by a first weighting factor before the mean is found, a z 
gradient which is situated in an overlap region of the object being multiplied by a larger first 
weighting factor than a z gradient that is not situated in an overlap region. An overlap region 
is a region of the intermediate image that has been reconstructed from measured values 
whose times of acquisition were situated in different periods. As mentioned at the beginning, 
10 motion artifacts occur to a more pronounced degree in these overlap regions, which means 
that a heavier weighting of the z gradients in the overlap regions will give an additionally 
improved reduction in the motion artifacts. For example, all the z gradients that are situated 
in an overlap region may be multiplied by 1 and all the z gradients that are not may be 
multiplied by 0. 

15 If there are known regions A of the object in which the object has a lower 

natural z gradient than in other regions B of the object, then a z gradient of the intermediate 
image that is situated in a region of the image corresponding to a region A of the object, may, 
in addition, be multiplied, before the mean is found, by a second weighting factor that is 
larger than a second weighting factor for a z gradient of the intermediate image that is 

20 situated in a region of the image corresponding to a region B of the object. In the extreme 
case, all the z gradients in regions of the image that correspond to region A of the object are 
multiplied by 1 and all the z gradients in regions of the image that correspond to region B of 
the object are multiplied by 0. This produces a reduction, by natural z gradients of the object, 
in the possible disruption of the z gradient caused by motion artifacts and hence to a further 

25 reduction of the motion artifacts. 

The z gradients may also be multiplied with a combination, and particularly a 
linear combination, of the first and second weighting factors before the mean is formed. What 
is more, two motion-artifact values, one obtained with the help of the z gradient and one with 
the help of the gradient in the direction of the angle of rotation of the radiation source, may 

30 be combined, in particular linearly, to give a single motion-artifact value. 

In addition to the motion-artifact value, a similarity value can be determined 
for each intermediate image. For this purpose, each intermediate image is compared with the 
intermediate images adjacent in time to give a similarity metric. The similarity metric may be 
any function that gives a value for similarity that is smaller the more similar two intermediate 
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10 

images are. The similarity metric may for example make use of correlations, and particularly 
the Pearson correlation, and differences in mutually corresponding image values in different 
intermediate images. A similar value a for two intermediate images may for example be 
obtained from the following equation: 

CT =^XK- F 2,,| where i = l,...,JV. (1) 

Here, N is the number of voxels in one of the intermediate images, Vyis the image value of 
the rth voxel in a first intermediate image and V^,; is the image value of the z'th voxel in a 
second intermediate image. The mean square deviation of mutually corresponding voxels 
may also be used as a similarity metric: 

CT =^ZV(^-^T where i = 1. - (2) 



Respective similarity values are determined for the first and last intermediate 
images in time and for all the other intermediate images two similarity values are determined, 
and in the latter cases a mean is formed of the two similarity values that are assigned to each 
intermediate image to give one similarity value, so that in the end precisely one similarity 

15 value is assigned to each intermediate image. 

If the similarity value of an intermediate image is relatively small, then this 
intermediate image differs to only a small degree from the intermediate images that adjoin it 
in time, which means that the corresponding object states too are not very different from one 
another. Hence the object hardly moved in the region of the phase point that is assigned to the 

20 intermediate image concerned. 

In step 106, an intermediate image is selected that has the lowest motion- 
artifact value. If, in addition, similarity values were determined in step 105, then for each 
intermediate image the similarity value can be combined with the motion-artifact value to 
give a combination value. The combination value may for example be a linear combination, 

25 and in particular the mean, of the motion-artifact value and the similarity value. What is then 
selected in step 106 is an intermediate image that has the lowest combination value. 

In step 107 the CT image is reconstructed, the measured values used being 
solely ones that were situated, in the particular period, in a phase of movement that was 
located at the phase point that corresponds to the intermediate image determined in step 106. 

30 The method according to the invention is not confined in this case to any particular method of 
reconstruction. The CT image may for example be reconstructed by means of a filtered back- 
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projection that is explained below in connection with Fig.4, or by means of known iterative 
methods. 

Before the reconstruction, each measured value may be weighted by 
multiplying it by a weighting factor, if this had not already been done when the 
5 corresponding intermediate image was being reconstructed in step 104. This weighting factor 
may be all the larger, the more the corresponding measured value is positioned in the center 
of the particular phase of movement. 

In contrast to the reconstruction of the intermediate images in step 104, the 
concluding reconstruction of the object or FOV in step 107 is performed with a high spatial 
10 resolution (e.g. 20 x 20 x 20 cm 3 represented by 5 12 3 voxels), in order to obtain an image of 
the highest possible quality. 

In what follows, a filtered back-projection by means of which the intermediate 
images and the concluding CT image can be reconstructed will now be described as an 
illustrative example (Fig. 4). 
15 For reconstruction, the measured values are regrouped in parallel in step 201 . 

What the parallel regrouping does is resort and re-interpolate the measured values as if they 
had been measured using a different radiation source (an extended radiation source that is 
positioned along part of a helix and is able to emit mutually parallel fans of rays in each case) 
and using a different detector (a plane, rectangular "virtual detector" containing the axis of 
20 rotation 14). 

This will be explained in more detail by reference to Fig. 5. Reference numeral 
1 7 in this Figure designates the helical trajectory from which the rays from the radiation 
source pass through the examination region. A fan-shaped bundle of rays 43, whose rays 
extend in a plane containing the axis of rotation 14, originates from the radiation source 
25 position So. The conical bundle of rays that is emitted by the radiation source at position So 
can be thought of as composed of a plurality of plane fans of rays that lie in planes parallel to 
the axis of rotation 14 and intersect at the radiation source position So. Fig. 5 shows only a 
single one of these fans of rays, namely fan of rays 43. 

Also shown in Fig. 5 are further fans of rays 41, 42 and 44, 45 that are parallel 
30 to the fan of rays 43 and lie in planes parallel to one another and to the axis of rotation 14. 
The associated radiation source positions 5-2, S-1 and Si, S2 are taken up by the radiation 
source S respectively before and after it reaches radiation source position S 0 . 

The fans of rays 41 to 45 form a group and define a bundle of rays 70 of a 
tent-like shape. A group of fans of rays is called a projection. For each projection, there is 
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now defined a rectangular virtual detector 160 that lies in a plane that contains the axis of 
rotation 14 and is oriented perpendicularly to the parallel fans of rays making up a projection. 
The points forming the corners of the virtual detector 160 are the points at which those rays 
that, starting from the outer positions of the radiation source, impinge on the section of the 
5 helix situated opposite pass through the said plane. For the bundle of rays 70 in Fig. 5, S. 2 and 
S2 are the outer radiation source positions. On the rectangular detector 1 60 are defined 
detector elements laid out to Cartesian coordinates, i.e. in rows and columns, to which the 
measured values are re-interpolated. 

Then, in step 202, the measured values assigned to the individual rays are 

10 multiplied by a weighting factor that is equal to the cosine of the included angle of the 

individual ray. The included angle of a ray is the angle that this ray makes with a plane that is 
oriented perpendicularly to the axis of rotation 14. When the said angle is small, its cosine is 
substantially equal to 1, and step 202 can therefore be dispensed with. 

In step 203, one-dimensional filtering with a transmission coefficient that rises 

15 in a ramp shape with spatial frequency is applied to the measured values. What are used for 
this purpose are values that succeed one another in a direction perpendicular to the axis of 
rotation 14, i.e. along a row of the detector 160. This filtering is performed along each row of 
the virtual detector for all the groups of fans of rays. 

In other embodiments, the parallel regrouping may be dispensed with. It is 

20 then known that the filtering has to be modified, because the detector unit is curved, e.g. in an 
arc, around the radiation source or the axis of rotation. 

In step 204, a voxel V(x) is determined within the FOV. Then, in step 205, a 
projection, i.e. a group of fans of rays, is selected that has not yet been used to reconstruct the 
voxel V(x) and whose time of acquisition was situated in one of the intervals of time defined 

25 above. If there is no ray in the projection that passes through the center of voxel V(x), it is 
determined at what point a ray that did pass through the center would impinge on the surface 
of the detector. The associated measured value is then calculated by interpolating from the 
measured values of adjacent rays. The measured value that can be assigned to that ray of the 
projection that passes through the voxel, or the corresponding measured value obtained by 

30 interpolation, is summed in step 206 to give the voxel V(x). In step 207, a check is made to 
see whether all the projections have been looked at. If they have not, the flow chart branches 
back to step 205. Otherwise, a check is made in step 208 as to whether all the voxels V(x) in 
the FOV have been dealt with. If they have not, the process continues with step 204. If on the 
other hand all the voxels V(x) in the FOV have been dealt with, then the absorption has been 
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determined in the whole of the FOV and the computer tomography method according to the 
invention is brought to an end at step 108 (see Fig. 2). 

Steps 201 to 203 can be dispensed with when reconstructing the concluding 
CT image if the measured values have already been treated in the appropriate way when the 
5 intermediate images were reconstructed in step 104. 

By the method according to the invention, an intermediate image having the 
lowest motion-artifact value or combination value can be determined for the entire object, i.e. 
the entire FOV. It is however also possible in accordance with the invention for the FOV to 
be divided into a plurality of sub- fields (sub-regions of the object), and for an intermediate 

10 image having a minimal motion-artifact value or combination value to be determined for 
each sub-field by steps 104 to 106, thus allowing an optimum phase point to be determined 
for each sub-field. In the concluding reconstruction in step 107, the measured values used for 
each sub-field are solely ones that were acquired while the object was in a phase of 
movement that was situated at the phase point determined for this sub-field. Treating the sub- 

15 fields differently in this way allows for the fact that the object may move differently in 
different ones of its regions, and enable motion artifacts to be further reduced in this way. 
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LIST OF REFERENCE NUMERALS 



h Distance between adjacent turns of a helical trajectory 

Ri, R 2 Points at which R peaks are situated 

S Radiation source 

5 S-2 . . . S2 Radiation source positions 

ct max Included angle 

q>\ ... (Pn Phase points 

1 Gantry 

2, 5 Motors 

10 3 Collimator arrangement 

4 Bundle of rays 

7 Control unit 

8 Electrocardiograph 
10 Reconstructing unit 

15 11 Monitor 

13 Examination region 

14 Axis of rotation 

16 Detector unit 

17 Helical trajectory 
20 21 Electrocardiogram 

23 R peak 

41 ... 45 Fans of rays 

70 Bundle of rays 

160 Virtual detector 
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CLAIMS: 



1 . A computer tomography method having the following steps: 

a) generation by a radiation source (S) of a bundle of rays (4) that passes though 
an object moving periodically, 

b) production of a relative movement between the radiation source (S) on the one 
5 hand and the object on the other hand, which relative movement comprises rotation about an 

axis of rotation (14), 

c) acquisition, with a detector unit (16) and during the relative movement, of 
measured values that depend on the intensity in the bundle of rays (4) on the far side of the 
object, 

1 0 d) sensing of a movement signal (21) dependent on the movement of the object 

with a movement-sensing means (8) and determination of cyclically repeated phases of 
movement with the help of the movement signal (21) sensed, 

e) reconstruction of a plurality of intermediate images of a region of an object, 
each intermediate image being reconstructed with measured values that were acquired while 

15 the object was in a different phase of movement, thus enabling a phase of movement to be 
assigned to each intermediate image, 

f) determination of the phase of movement in which there was least movement of 
the object in the region, by determining the intermediate image having the fewest motion 
artifacts in the region, 

20 g) reconstruction of a computer tomographic image of the region from measured 

values that were acquired while the object was in the phase of movement in which there was 
least movement of the object in said region, the reconstruction parameters that are used in 
this case differing from the reconstructions parameters used to reconstruct the intermediate 
images. 

25 

2. A computer tomography method as claimed in claim 1 , characterized in that 
the intermediate images in step e) are reconstructed with a lower spatial resolution than the 
computer tomographic image to be reconstructed in step g). 
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3. A computer tomography method as claimed in claim 1 , characterized in that 

the region of the object that is to be analyzed (the FOV) is divided into a plurality of sub- 
regions and in that steps e) to g) are performed for each sub-region. 

5 4. A computer tomography method as claimed in claim 1 , characterized in that, 

with the help of a motion-artifact metric, there is determined for each intermediate image a 
motion-artifact value by applying the motion-artifact metric solely to measured values from 
the particular intermediate image, and in that the intermediate image having the lowest 
motion-artifact value is determined to be the intermediate image having the fewest motion 
10 artifacts. 



5. A computer tomography method as claimed in claim 4, characterized in that 
the motion-artifact value of an intermediate image is the mean of gradients of image values in 
the intermediate image in the direction of an axis of rotation. 

15 

6. A computer tomography method as claimed in claim 5, characterized in that 
the gradients are weighted before a mean thereof is formed, in which case a gradient that is 
situated in an overlap region of the object, through which region rays having acquisition 
times situated in different periods pass, is given a higher weight than a gradient that is not 

20 situated in an overlap region. 



7. A computer tomograph for carrying out the method claimed in claim 1 , having 

a radiation source (S) for generating a bundle of rays (4) that passes 
through an object moving in a cycle, 
25 - a drive arrangement (2, 5) for producing a relative movement 

between the radiation source (S) on the one hand and the object on the other hand, which 
relative movement comprises a rotation about an axis of rotation (14), 

a detector unit ( 1 6) for acquiring, during the relative movement, 
measured values that depend on the intensity in the bundle of rays (4) on the far side of the 
30 object, 

a movement-sensing means (8), in particular an electrocardiograph 
(8), for the sensing of a movement signal (21) dependent on the movement of the object with 
a movement-sensing means (8), 

a reconstructing unit ( 1 0) for reconstructing a computer tomographic 
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image of the object from the measured values, 

a control unit (7) for controlling the radiation source (S), the drive 
arrangement (2, 5), the detector unit (16), the movement-sensing means (8) and the 
reconstructing unit (10) in the following steps: 
5 a) generation by the radiation source (S) of a bundle of rays (4) that 

passes though an object that moves in a periodically, 

b) production of a relative movement between the radiation source (S) 
on the one hand and the object on the other hand, which relative movement comprises 
rotation about an axis of rotation (14), 
1 0 c) acquisition, with the detector unit ( 1 6) and during the relative 

movement, of measured values that depend on the intensity in the bundle of rays (4) on the 
far side of the object, 

d) sensing of a movement signal (2 1 ) dependent on the movement of the 
object with the movement-sensing means (8) and determination of periodically repeated 

1 5 phases of movement with the help of the movement signal (21) sensed, 

e) reconstruction of a plurality of intermediate images of a region of the 
object, each intermediate image being reconstructed with measured values that were acquired 
while the object was in a different phase of movement, thus enabling a phase of movement to 
be assigned to each intermediate image, 

20 f) determination of the phase of movement in which there was least 

movement of the object in the region, by determining the intermediate image having the 

fewest motion artifacts in the region, 

g) reconstruction of a computer tomographic image of the region of the 

object from measured values that were acquired while the object was in the phase of 
25 movement in which there was least movement of the object in said region, the reconstruction 

parameters that are used in this case differing from the reconstruction parameters used to 

reconstruct the intermediate images. 

8. A computer program for a control unit (7) for controlling a radiation source 

30 (S), a drive arrangement (2, 5), a detector unit (16), and a reconstructing unit (10) of a 
computer tomograph for carrying out the method claimed in claim 1 . 
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